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We present a comparative analysis of the velocity dependence of atomic-scale friction for the 
Tomlinson model, at zero and finite temperatures, in ID and 2D, and for different values of the 
damping. Combining analytical arguments with numerical simulations, we show that an appreciable 
velocity dependence of the kinetic friction force Ff r i C , for small scanning velocities v 3 (from 1 nm/s to 
2 fim/s), is inherent in the Tomlinson model. In the absence of thermal fluctuations in the stick-slip 
regime, it has the form of a power-law, Ff r i C — Fo oc vf with j3 — 2/3, irrespective of dimensionality 
and value of the damping. Since thermal fluctuations enhance the velocity dependence of friction, we 
provide guidelines to establish when thermal effects are important and to which extent the surface 
corrugation affects the velocity dependence. 
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I. INTRODUCTION 



Although friction is a common phenomenon in every- 
day experience, the fundamental mechanisms governing 
friction at the atomic level are still under discussion. 
For macroscopic contacts the friction force is found to 
be independent of the sliding velocity, but no consen- 
sus has been reached on the velocity dependence at the 
nanometer scale. A very powerful technique for mea- 
suring atomic-scale friction is provided by Atomic Force 
Microscopy (AFM) [1,2]. Since scanning velocities ac- 
cessible by AFM are very small, typically from nm/s to 
few //m/s, it is relevant to study friction dynamics in 
this regime. Velocity dependence of friction is relevant 
both for applications and from a fundamental point of 
view, and has been discussed in several AFM [3-11] and 
Quartz Crystal Microbalance [12] experimental studies 
as well as theoretical works [9-11,13-19]. Depending on 
the investigated systems and on the experimental con- 
ditions, different and somewhat contradictory results for 
the velocity dependence have been found. In the original 
experiments of Mate et al. [3] the authors state that the 
frictional forces of a tungsten tip on graphite show lit- 
tle dependence on velocity for scanning velocities up to 
400 nm/s. A similar behavior up to velocities of several 
/xm/s has been reported also in the work of Zworner et 
al. [10], where friction on different carbon structures has 
been studied. The authors of Ref. [10] claim that a ID 
Tomlinson model at T — can reproduce a velocity in- 
dependent friction force for scanning velocities up to ~ 1 
/xm/s, while giving a linear increase of friction for higher 



velocities. At variance with the ID case, in the 2D ver- 
sion of the Tomlinson model at T = 0, which has been 
recently analyzed by Prioli et al. [11], a smooth increase 
of friction for velocities lower than ~ 300 nm/s has been 
found. In view of the results of Zworner et al. for the 
ID case, the authors argue that this effect should be pe- 
culiar of the 2D model, due to the non-linear coupling 
between the two degrees of freedom in the system. The 
role of damping has not been addressed in Refs. [10,11]. 
In the underdamped regime, the velocity dependence can 
be quite complex, especially at intermediate-large veloci- 
ties, where the system displays bifurcations, chaotic mo- 
tion, resonances and hysteresis [14]. In the overdamped 
regime, Robbins and Miiser [20] suggest velocity inde- 
pendent friction. 

An increase of the friction force has been observed for 
small velocities also in Refs. [6,7,9] and it has been at- 
tributed to thermally activated processes [6,7,9,19]. By 
means of a simple thermal activation probabilistic anal- 
ysis in ID, Gnecco et al. [9] have obtained a logarithmic 
increase of friction with scanning velocity which fits their 
experimental data quite well. A similar dependence had 
been obtained using a simple stress-modified thermally- 
activated Eyring model [6]. In a recent work, Sang et 
al. [19] have corrected this logarithmic relation at not too 
small velocities: they propose a | lnw s | 2 / 3 dependence of 
the friction force, where v s is the scanning velocity. How- 
ever, recent experiments showing an increase of friction 
with velocity [11] do not display the logarithmic behav- 
ior related to thermal activation, but rather suggest an 
athermal power-law v@ behavior, as found in related sys- 
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terns, such as charge density waves [21] and in boundary 
lubrication [22]. 

In view of the contradictory results presented above, 
here we reexamine this issue for Tomlinson-like models 
in ID and 2D, for different values of the damping, and 
both with and without thermal effects. In particular, we 
focus on the importance of the athermal contribution to 
the velocity dependence of friction, which is intrinsically 
present in the Tomlinson model. We show by means of a 
combined analytical and numerical analysis that the ex- 
ponent (3 is independent of the spatial dimension and of 
the damping. Then we discuss the role of thermal fluc- 
tuations, establishing guiding rules to understand where 
thermal effects become dominant. 

In Sec. II we illustrate the model studied and the nu- 
merical techniques. In Sec. Ill we discuss the results 
for the athermal velocity dependence of friction and in 
Sec. IV we include thermal fluctuations. Finally, we 
present some concluding remarks in Sec. V. 



II. MODEL 

The Tomlinson model [23] has been successfully used 
to describe the motion of a tip and to model the scan 
process in AFM [24-27]. In particular, this model can 
reproduce the stick-slip motion observed in experiments 
and can be used to study frictional dynamics. Here we 
consider the ID Tomlinson model and its extension in 
2D at T — and T ^ 0. A cantilever tip of mass m in- 
teracts with the surface via a periodic potential V ts and 
is attached by a spring of elastic constant k x to a support 
moving at constant velocity v s along the x direction. For 
the ID case we choose V ts of the form 



V ts {x) = V [l - cos(2wx/a x )}, 



(1) 



where a x is the lattice constant of the substrate. The 
clastic interaction between the tip and the support is 



Vei(x) = ^k x (x-x s 
where the support position x s is 



v s t. 



(2) 



(3) 



It is assumed that the tip is a point-like object, represent- 
ing the average over many atoms of the real tip-surface 
contact. Energy dissipation in this model is introduced 
by adding a damping term proportional to the tip veloc- 
ity in the equation of motion. Thermal fluctuations are 
taken into account by a stochastic force, in the frame- 
work of the Langevin approach. Thus, the equation of 
motion in ID becomes 



mx 



2itVq . f 2irx\ 
+ rnqx H sin I J + k x {x - v s t) = fit), (4) 



with the random force f(t) satisfying the conditions 

< f(t) > = and < /(t)/(0) >= 2mr}k B T5{t), where 

< • > indicates an ensemble average, r\ is the damping 
parameter and fee is the Boltzmann's constant [28]. The 
static friction force in this model is simply given by the 
force needed to overcome the potential barrier: 



F, 



2ttV 



static 



(5) 



Now we discuss the behavior of the ID model at T = 0, 
i.e. when f(t) = in Eq. (4). In this situation the 
solution of Eq. (4) for T = is periodic, with period 
na x /v s [14]: 

x(t + na x /v s ) — x(t) + na x for integer n. (6) 

Usually n = 1 for not too small r\. 

Elastic instabilities leading to nonadiabatic jumps be- 
tween metastable states occur for soft cantilever spring 
constants, in particular when [24,27] 



k x <C 



d 2 Vu 



dx 2 



i.e. Vb 



4tt 2 Vc 



k x a x 



> 1, 



(7) 



where x m = na x denotes the position of the minima of 
Vts- In this case stick-slip motion, often observed in AFM 
experiments, is expected and the kinetic friction force is 
finite in the limit v s — > 0. Conversely, for Vo < 1, uni- 
form sliding occurs and energy dissipation comes only 
from the viscous term mriv s , which vanishes for v s — > 0. 
Notice that the kinetic friction force for v s — > is not 
equal to the static friction force F sta tic, since it results 
from dynamical effects and not by the interaction poten- 
tial Vts. The kinetic friction force Ff r i C is defined as the 
mean value of the lateral force F x = k x (v s t — x) over 
time [10,14,27]. By assuming a periodic motion of the 
type of Eq. (6) , Ff r i C can be written as 



Ffric — ^ F x 



Vs f 

na x J 



na x jv s 



F x dt. 



(8) 



It is easy to show that the definition Eq. (8) is equivalent 
to calculating the friction force from the energy dissipa- 
tion AW in one period 



AW 



— mr) 

Jo 



na x /v s 



x 2 dt. 



The friction force is given by 

L fne 



na 7 



(9) 



(10) 



Here we extend the model to deal with the motion at 
zero and finite temperature on a 2D lattice, as done in 
Refs. [11,27] for T = 0. The tip-surface interaction is 



/2ttx\ [2-Ky\ 
V ts (x,y) = Uocos y—J cos [j^J ' 



(11) 
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where a x and a y are the lattice parameters in the x and 
y directions respectively. When a y = \/3a x the substrate 
has the symmetry of a hexagonal closed-packed lattice. 
The clastic interaction is 



V e i(x,y) = -k x (x - v s tf + -k y (y - y s 



(12) 



where k y denotes the spring constant in the y direction 
and y s = constant represents the scanning line of the 
support. The equations of motion can be written in 2D 
as 



rax + mr}x - V sin cos (^p) + k x (x - v s t) 

my + m-qy - V cos (^fj sin J + k y (y - y s ) 



/*(*) 

hit) 



may affect the dynamical be havior of the system. Usually 
a critical damping, 77 = 2y/k x /m [27], is assumed. Here 
we study the problem for different values of 77, in the un- 
derdamped, overdamped and critically damped regime. 
For each fixed scanning velocity v s , we compute the fric- 
tion force Ff r i c , averaging over many stick-slip periods 
(usually 10 at T = and 100 at T ^ 0), according to 
Eqs. (8) and (14). The behavior of Ff r i C as a function of 
v s in ID is shown for critical damping in Fig. 1(a) on a 
linear scale and in Fig. 1(b) on the most commonly used 
log- log scale [10]. Notice that the log-log scale hides the 
velocity dependence for small velocities (v s < 1.5 /im/s), 
where the friction force varies by more than 10%. 



(13) 



where f x and f y are independent stochastic forces satis- 
fying the same properties as / in Eq. (4) . In this case we 
also have a component of the lateral force along y, i.e. 
F y = k y (y s — y). The definition of the friction force in 
Eq. (8) can be generalized in 2D as 



?fric = ^/<F x > 2 + <F y >' 



(14) 



We have solved the non-linear equations (4) and (13) us- 
ing a Runge-Kutta 4 algorithm with initial conditions 



x(Q) = 0, i(0)=0, 2/(0) = 0, y(0) 



0. 

(15) 



and for different values of the scanning velocity v s and of 
the scanning line y s . 



III. ATHERMAL VELOCITY DEPENDENCE OF 
FRICTION 

At T = the dynamics can be described by the 
equations of motion (4) and (13) without the stochas- 
tic forces. We choose values of the parameters which are 
typical of AFM experiments: m — 10~ 10 kg, k x = 10 
N/m [7,27,29], a x = 0.316 nm (in 2D we set a y = 0.548 
nm, corresponding to the hexagonal-packed structure of 



MoS 2 (001) [27], and k x 



giving a resonance fre- 



quency \Jk x /m of the order of 10 5 Hz, which is charac- 
teristic of AFM experiments. In principle, the corruga- 
tion Vq of the tip-surface potential depends on the load- 
ing force, which is not considered in ID and 2D models: 
typically Vq ranges from 0.2 eV to 2 eV, as found in dif- 
ferent studies [30,31]. Thus we take V = 1 eV. These 
values of the parameters give Vq = 7, yielding stick-slip 
motion (Vo > 1) and allowing us to compare directly our 
results with those of Zworner et al. in ID [10]. The time 
step used in the calculations is ~ 0.1 ns, a value which is 
needed to account for the fast oscillations in the under- 
damped regime. The choice of 77 is quite delicate and it 
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FIG. 1. Frictional force Ff ric as a function of sliding veloc- 
ity v 3 in the ID Tomlinson model, plotted on a linear (a) and 
on a log-log scale (b) for Vq = 1 eV, m = 10~ 10 kg, k x = 10 
N/m, a x — 0.316 nm (Vq = 7) and r\ — 2y // k x /m ~ 6.3 ■ 
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s The increase of Ff r i c for small velocities is hidden using 
a log-log scale. The dotted line in (a) is a power-law fit to the 
data of the form Ff ric — F a oc vV A for v s < 2 /_tm/s. 

The data in Fig. 1 (a) can be fitted quite accurately by 
a power law of the form 



F 



fric 



Fn 



(16) 



with P ~ 2/3 and c a constant depending on the param- 
eters of the model and on the space dimension. 

To our knowledge the athermal velocity dependence of 
atomistic dry friction has been scarcely investigated up 
to now: it has been studied in the limit of large veloci- 
ties [14] and in the context of boundary lubrication [22] . 
Here we discuss the velocity dependence of dry friction for 
small scanning velocities, in the stick-slip regime, which 
is described by Eq. (16). In this case, the value of the 
exponent (3 can be calculated analytically for the Tom- 
linson model, yielding (3 — 2/3, as we will show below. 
The same kind of behavior has been found in the field 
of elastic manifolds, for the dynamics of charge density 
waves driven by an electric field [21] and for the motion 
of a contact line on a heterogeneous surface [32,33]. This 
law characterizes the athermal motion of strongly pinned 



3 



systems (Vb > 1 in our terminology), moving at constant 
velocity. 

Considering for simplicity the ID case and following 
Rcf. [21], we look for a solution x(t) of Eq. (4) in the 
athermal case (/(t) = 0) of the form 



Since energy dissipation takes place mostly near the fast 
slip events, we focus on the behavior of Eq. (25) near the 
slip point x , determined by 



x(t)=x A (t) + 6(t), 



(17) 



dVtot , , , . 2tt . 
— — = k x (x - x s ) H Vb sin 

ax a x 



/27rar\ _ 
\ a x ) 



where xa is the adiabatic solution of Eq. (4), i.e. the so- 
lution for v s — > 0, and is a perturbation. The adiabatic 
solution satisfies Eq. (4) neglecting the first (inertial) and 
second (damping) term: 



d 2 Vto 
dx 2 



— k x 4- 



Vq cos 



2ttx\ 
a x J 



= 0. 



(26a) 
(26b) 



where Vj t = V< s + V e i is the total potential energy. From 
Eq. (26b) the position x of the tip right before a slip 
event is 



k x (x A - v s t) = sin 

a x \ a x J 

From Eq. (8) it follows that 



(18) 



x Q = arccos(Vb). 

Z7T 



(27) 



-F/ric = < fc x (w s i - xa - 0) >= 

k x < (v s t - x A ) > -k x < 9 >= F - k x < 9 >, (19) 

having defined F = Ff ric (v s — > 0). Thus, the final goal 
is to work out the dependence of 



Eq. (26a) gives the position x^ of the support at the 
slip point: 



T (0) _ ^£ 
^ ~2^ 



1 + arccos 



1 



na x /v s 

<9>=— I e(t)dt 



v. 

na x jo 



(20) 



Near the slip point we can set 

x A (t) = x +£ A (t) 

with 



on v s . First we notice that for Vb ^ 1 the inertial term 
mx can be neglected with respect to the damping term 
mrji near a slip event. This can be straightforwardly seen 
in the adiabatic limit. In fact, differentiating Eq. (18) 
with respect to time we obtain 



\u\ « 



2tt' 



Using Eqs. (7) and (26b) and the relations 
'2ttx 



k x x A k x v s — 
giving for x A and x A 

Z A = X A = — 



Vq COS 



k x v s 



( 2irx A \ . 

xa, 

\ a x J 



sin 



(21) 



/ 2ttxa \ 
\ a x J 



sin 



+ 



2tt\ 
a x J 



cos 



and 



x A = z A 



dz A 

dXA 



ZA = 



fe) Vbcos(^) 
{k x v s f (fr) 3 ^osin(^) 



(22) 



(2ttx a \ 

cos ~ 

\ a x J 

Eq. (25) becomes 



cos 



2irx 



sm 



27TXr 



2ttx \ 



(28) 

(29) 
(30) 

u 
u 



1 e (vj - 1)64 



<C 7/. 



(31) 



^ + (!r) 2 ^cos(2^) 



Since Eq. (30) holds we can neglect the second term with 
respect to the first, obtaining 



respectively. Then the condition 
\x A \ < v\xa\ 

becomes 

k x v s Vo(%y 



sm 



^ + (fr) 2 ^cos(2^) 



(23) 
(24) 

(25) 



1/2 



\u\ » 



v s a x 



2^^Vf- 



(32) 



1 



Eq. (32) is easily fulfilled for large Vb (or large 77) and/or 
small v s . For example, with our choice of parameters, 
yielding Vb ~ 7, and 77 ~ 6 • 10 5 s _1 , conditions (32) is 
valid for velocities up to v s <~ /um/s. Having now demon- 
strated that we can neglect the inertial term, we can 
expand the equation of motion (without the term mx) 
near xq: 



4 



mvi = k x v s St + i ( — 
2 \a x 



Vq sin 



/ 2irx \ 
\ a x J 



where 



and 



£ = x - x 



St = t — to, 



(33) 



(34) 



(35) 



to being the time at which the slip takes place. Following 
Rcf. [21], with the change of variables 



where C = ^ 



St = Cv7 l ' Z T 



_Vi_ sin (W) 

2mr) V a x I 



-1/3 



form of a Riccati equation: 

dx 2 
Tt =X - 



m7y 



-T. 



(36a) 
(36b) 

, Eq. (33) takes the 
(37) 



It can be shown [21] that the major contribution to the in- 
tegral (20) comes from a time St = St s = t\ — to such that 

— 1 /3 

8t s oc v s When t ~ t\ the solution x{ T ) OI the Ric- 
cati equation has a divergence of the form x( T ) ~ T L T ■ 
Note that 5i s is the slip time, i.e. the time it takes for 
the tip to go from the metastable position x = x to the 
next metastable position x = x\. For the adiabatic so- 
lution the slip occurs instantaneously, while St s is finite 
for finite v s and this is responsible for the velocity depen- 
dent correction of the friction force. In fact, when t <~ t\ 
£ ~ x\ — xq is of order 1 (e.g. independent of v s ), and 
9 = x — xa = £, — S,a is of order 1 as well. Thus 



< 



6(t)dt (x v 2 J 3 , 



(38) 
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FIG. 2. Tip position as a function of support position in 
the ID Tomlinson model for different values of the scanning 
velocity (from left to right v s = 1.5 nm/s, 15 nm/s, 300 nm/s, 
750 nm/s, 1.5 /im/s), r\ — 2yjk x /m and Vo = 7. The inset is 
a blow up of the region around the first slip event. 

Furthermore the slip position tends to move right- 
wards for increasing v s This means that the integral of 
F x = k x (x s — x) over one period 



~" fric 



i r 

na x J 



F x dx s 



k x {na x f 



na x 



k x f 
na x 7 



(39) 



*x J t 



increases with increasing v s , since the second term on the 
right side of Eq. (39) decreases. Fig. 3 shows the slip time 
St s as a function of v s , as measured from the numerical 
solution of the equation of motion. 



0.20 , 



0.16 



0.12 



0.08 



which proves that the exponent /3 appearing in Eq. (16) 
is (3 — 2/3. This shows that the dependence of friction 
on velocity is a dynamical effect which is due to the fi- 
nite (although small) scanning velocity, as it can be seen 
in Fig. 2, where the tip position X cLS cL function of the 
support position x s is plotted. The important feature is 
that the slip events are not instantaneous, as highlighted 
in the inset of Fig. 2, showing a finite slip time which 
depends on v s . Only if the slip events were really instan- 
taneous a velocity independent friction force would natu- 
rally follow from the definition Eq. (8), giving Ff ric = F . 
Therefore, the source of athermal velocity dependence of 
friction is the non adiabaticity of the motion of the tip 
for finite u,. 
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FIG. 3. Slip time as a function of scanning velocity in the 
ID Tomlinson model for critical damping and Vo = 7. The 
points connected by the solid line are obtained by numerical 
simulations, while the dotted line is a power-law fit to the 
data of the form St s oc u^ 1 ^ 3 - 

The behavior of St s is in very good agreement with the 
scaling relation 



St s oc v- 1 '*, 



(40) 



•5 



which is the law expected from the discussion following 
Eq. (37). 



A. Effect of damping 

The effect of the damping parameter on the velocity 
dependence of friction has not been investigated so far 
in the literature, because the typical choice is to assume 
critical damping in order to damp out the fast oscillations 
of the tip after the slip events and to avoid jumps of the 
tip of more than one lattice parameter. Nevertheless, it 
would be desirable to know the dynamical behavior of 
the tip for a range of values of 77, since experimental situ- 
ations do not always meet the condition of critical damp- 
ing. The behavior of Ff r i C vs. v s , for values of 77 ranging 
from strongly underdamped to strongly overdamped, is 
reported in Fig. 4. 




3=2 
ti=1.5 



Tl=0.4 - 



dent friction contrary to the qualitative expectation of 
Ref. [20] . The authors of Ref . [20] argue that in the over- 
damped regime the peak velocity of the tip, correspond- 
ing to a slip event, is a constant equal to 2ttVq / (mrja x ) . 
This would imply that the amount of energy dissipated, 
which is proportional to the tip velocity according to 
Eq. (9), should not depend on v s . On the contrary, we 
find appreciable dependence also in this case. As it can 
be seen from Fig. 5, the peak velocity of the tip is not a 
constant, but increases appreciably by increasing v s . 

2.5, , , , , , 
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FIG. 5. Tip velocity as a function of support position in the 
ID Tomlinson model for different scanning velocities (from 
left to right v 3 = 1.5 nm/s, 15 nm/s, 300 nm/s, 750 nm/s) in 
the overdamped case (77 = 100y / kjm) and for V = 7. The 
horizontal line is the value 2nVo / (77177 a x ). 
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FIG. 4. Frictional force Ff ric as a function of sliding ve- 
locity v 3 in the ID Tomlinson model for Vo — 7 and differ- 
ent v alues of the damping parameter: from bottom to top 
v/(y/kx/m) = 0.4, 1.5, 2, 10, 100. The dotted lines are fit to 
the numerical data of the form Ff ric — Fa oc v% , with /3 = 2/3. 
In the most underdamped case (lower line) the friction force 
is lower because the tip performs jumps of two lattice param- 
eters. 

All curves start from the same value F , except for 
very low 77 (see discussion below), and can be fitted by 
Eq. (16) with the same value of /? = 2/3, suggesting that 
the functional form of the velocity dependence of friction 
is robust with respect to the strength of the damping. 
The velocity range of validity of Eq. (16) decreases for 
large 77, because the viscous regime (Ff ric ~ mriv s ) sets 
in for smaller values of v 8 (for example the data in Fig. 4 
are fitted up to v s — 1.2 /zm/s for 77 = 2^k x /m and up 
to v s — 0.3 jum/s for 77 = lQQ\Jk x /m). As expected, at a 
fixed value of v s > 0, Ff ric increases with 77, since energy 
dissipation increases by increasing the damping (see also 
Eq. (9)). Moreover the value of c in Eq. (16) is larger for 
larger 77, reflecting the fact that the variation of Ff r i C is 
more pronounced for the highest values of 77. 

Note that for high damping we find velocity depen- 



The lower curve in Fig. 4, corresponding to the highly 
underdamped value 77 = 0.4, is characterized by a much 
lower friction force, because the tip in this case can per- 
form jumps with periodicity of two lattice parameters 
(i.e. n = 2 in Eq. (6)). This makes the lateral force drop 
to lower values after a slip event with respect to the criti- 
cally damped situation, as shown in Fig. 6, resulting in a 
smaller F . Notice that in Fig. 6 we also plot the so called 
"mechanistic Tomlinson loop", i.e. F x = s i n (*™\ 

as a function of X j clS obtained from Eq. (26a). The slip 
events correspond to transitions between stable branches 
of this loop. 



B. Role of dimensionality 



As already mentioned in the introduction, this problem 
was recently studied in Ref. [11] using a 2D Tomlinson 
model, where a velocity dependent friction force was ob- 
served even for scanning velocities less than 300 nm/s. 
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FIG. 6. Lateral force as a function of tip position for 
two val ues of the damping parameter: critically da mped, 
rj = 2W k x /m (solid line) and underdamped, rj — 0.4^/ k x /m 
(dashed line). The reduced corrugation is Vo = 7 and the 
scanning velocity v s — 300 nm/s. Notice the presence of 
jumps with periodicity 2a x in the underdamped case. The 
upper and lower horizontal lines indicate the friction force 
for T) — 2^Jk x /m (Ff ric — 2.33 nN) and r\ = 0.4^/ k x /rn 
(Ff r i C — 1.01 nN) respectively. The dotted line represents 
F x = 2lla sin (^), as obtained from Eq. (26a). 




1.5 2.0 2.5 

x [nm] 

FIG. 8. Trajectory of the tip in the 2D Tomlinson model 
for critical damping, Vo = 7 and v s = 7.5 nm/s. The circles 
connected by the solid line indicate the positions of the tip in 
the xy plane during the dynamics. The regions where the dis- 
tribution of points is denser are the sticking domains, where 
the tip stays predominantly for most of the time. Note that 
the tip slips from one sticking domain to the other following 
a zig-zag pattern around the scanning direction (indicated by 
the dashed line, y s — 0.137 nm). 



Since for ID motion no velocity dependence had been 
previously found in Ref. [10], the authors attributed this 
dependence to the coupling between the two degrees of 
freedom of the system. Our results for the ID Tomlin- 
son model already give a dependence on velocity, and it 
is interesting to look at the effect of an extra dimension 
on this dependence. Indeed, as it can be seen in Fig. 7, 
the behavior of Ff r i C vs. v s in 2D for different values of 
the scanning direction y s is very similar to that in ID. 




1 ' 2 1 2 3 4 5 6 7 8 
v s [nm/s] 

FIG. 7. Friction force as a function of scanning velocity 
in ID (upper curve) and 2D Tomlinson model, for critical 
damping, Vo = 7 and different values of y 3 (from bottom to 
top y s = 0.274 nm, 0.137 nm, 0.1 nm and 0.05 nm). 

Thus, in spite of the 2D character of the tip motion, 
clearly visible in Fig. 8, no dramatic effect of the dimen- 
sionality on the velocity dependence of friction can be 
noticed. 



This result is actually not surprising, because the Tom- 
linson model is a mean-field model and the functional 
form of constituent relations, such as Ff r i C (v s ) should not 
change with dimensionality. Thus Eq. (16) is expected 
to hold also in 2D, with the same exponent j3 = 2/3. 
The values of the parameters Fq and c can be different 
in ID and 2D. Specifically Fo is always lower in 2D. In 
fact, in ID the tip is necessarily moved along an atom 
row , while in 2D, depending on the scanning line y s , 
the motion of the tip can occur also between atom rows. 
For the hexagonal lattice we have chosen, the interac- 
tion between the tip and the surface is the weakest when 
y s = a y /A (bottom curve of Fig. 7), while it reaches its 
maximum value for y s = 0, which coincides with the ID 
case (upper curve of Fig. 7) . Since the corrugation of the 
tip-surface interaction is directly related to the friction 
force [31], different scanning lines result in different val- 
ues of friction. This feature allows for example to obtain 
2D surface maps in AFM experiments (see for example 
Ref. [2]). We notice that the absolute variation of Ff ric 
with velocity in the lowest curves of Fig. 7 is more pro- 
nounced, thus supporting to a certain extent the claim of 
Ref. [11]. But it is important that this variation is only 
due to the different values of the prefactor c in Eq. (16) 
and not to a change of the exponent (3. Therefore, we can 
argue that no qualitative differences arise in the velocity 
dependence of friction in the 2D case and that the com- 
mon mechanism which produces the observed behavior 
at T = can be ascribed to the delayed athermal motion 
of the tip with respect to the support. 
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IV. EFFECT OF THERMAL FLUCTUATIONS 

At finite temperature we integrate numerically the full 
equations of motion Eqs. (4) and (13). Due the presence 
of the stochastic forces, the motion of the tip is quite 
noisy and averages over long trajectories (containing up 
to 100 periods) have to be considered in order to have a 
reliable value of the friction force. A typical behavior of 
the lateral force in ID for different velocities and critical 
damping at T = 300 K is displayed in Fig. 9. 




x [nm] 

FIG. 9. Lateral force as a function of tip position in the ID 
Tomlinson model for critical damping, T — 300 K and Vo = 7, 
for different scanning velocities (non-solid lines from bottom 
to top v s = 1.5 nm/s, 15 nm/s, 300 nm/s, 750 nm/s). The 
solid line represents F x = 2 ^ v ° sin (^^) , as obtained from 
Eq. (26a) (see also Fig. 6) . The inset shows a blow up of the 
region around a slip event. 

The height of the maximum for a fixed v s is not con- 
stant and the effect of the scanning velocity on the posi- 
tion of the slip is rather pronounced even for small v s . In 
fact, theoretical investigations based on simple analyt- 
ical approaches in II? [9,19] and numerical simulations 
of the ID Tomlinson model at T ^ [19] have shown 
that temperature is effective in overcoming the energy 
barriers AE, activating jumps of the tip between min- 
ima of the total potential energy, for temperatures such 
that AE ~ ksT. The thermal activation gives rise to a 
linear logarithmic dependence of friction on velocity for 
very small scanning velocities [9] : 

Ffric -F c k ln(v s ). (41) 

For a larger range of v s the following functional form has 
been proposed [19]: 

F fric -F c ot |ln(v s )| 2 / 3 . (42) 

The constant value F c is the lateral force corresponding 
to a slip event at T = 0. Eq. (42) is obtained by assuming 
Vo > 1 and Vq >• ksT . As it is shown in Fig. 10, where 
we compare Ff r i C vs. v s for T = and T = 300 K, the 
main source of velocity dependence of friction is due to 
thermal fluctuations in the system. 



3.0 r 
2.8 - 



2.6 - 




1 6 001 0.01 0.1 1 10 

v s [nm/s] 

FIG. 10. Velocity dependence of friction force in the ID 
Tomlinson model at T = (upper curve) and T = 300 K 
(lower curve) for critical damping and Vo = 7. The solid line 
is a fit of the data for T = 300 K, using Eq. (42) in the small 
velocity regime (v s < 2 /im/s). 

The data for T = 300 K can be fitted by a logarithmic 
behavior with exponent which is very close to the value 
2/3ofEq. (42). To our knowledge theoretical approaches 
of velocity dependence of friction at finite temperature 
have been restricted to ID models. Here we report re- 
sults of numerical simulations also for the 2D Tomlinson 
model, using the same parameters as for the model at 
T = 0. Not surprisingly, Fig. 11 shows that the velocity 
dependence of friction is very similar in ID and 2D, as 
we have found for T = 0. We can use Eq. (42) to fit the 
data of the 2D model as well. 

3.0 , 




0.001 0.01 0.1 1 10 

v s [nm/s] 

FIG. 11. Velocity dependence of friction force in the ID 
(upper curve) and 2D (lower curve) Tomlinson model for 
T — 300 K, critical damping and Vo = 7. The solid lines 
are fits to the data using Eq. (42) in the small velocity regime 
(v s < 2 /im/s). 

In fact, as we have discussed in Sec. IIIB, the mean 
field character of the Tomlinson model, preserves the 
same form of the velocity dependence of energy dissi- 
pation. 

The different behavior of the friction force with scan- 
ning velocity at T ^ is due to the activated motion of 
the tip, which lowers the friction force with respect to 
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the athermal situation. This can be easily understood 
from a sketch of the evolution of the total potential V to t 
during the scanning, which is given in Fig. 12. 




0.15 0.20 0.25 0.30 
x [nm] 

FIG. 12. Total potential energy Vtot as a function of tip 
position x for three values of the cantilever position x 3 (from 
bottom to top x s = 0.287 nm, 0.382 nm, 0.413 nm). The hor- 
izontal lines indicate the values of the minimum (V m in) and 
the maximum (Vmax) of the potential for each curve. The 
potential barrier is AE = Vmax — Vmin- The upper curve 
corresponds to AE — 0, while the middle curve to the case 
where AE ~ k B T. 



While at T = a slip event can occur only when the 
energy barrier AE (i.e. the difference between the max- 
imum and the minimum of Vtot) vanishes, thermal fluc- 
tuations can activate jumps of the tip from a metastable 
minimum to the next even for finite AE, when the can- 
tilever has reached a position which is smaller than the 
one needed for a slip at T = 0: specifically thermal ef- 
fects start to be significant as soon as AE is few times 
ksT. This has the effect to lower the energy dissipated 
in a jump, and thus the friction force. The energy barrier 
is given by 



max (*)) - V to t(Xmin(t)), 



(43) 



where x m i n and x max are respectively the positions of a 
metastable minimum and maximum of Vtot ■ 

Fig. 13 compares the velocity dependence of the fric- 
tion force for three values of Vq in the stick-slip regime 
(Vo = 0.28 eV, 0.57 eV and 1 eV), with k x = 10 N/m 
(yielding Vo = 2, 4 and 7 respectively), both for T = 
and T — 300 K. At the smallest scanning velocity con- 
sidered, in going from T = to T = 300 K, Ff ric de- 
creases only by a factor 1.2 for Vq = 7, but by a fac- 
tor 15 for Vo = 2. Indeed, by increasing Vo, the friction 
force Ff r ic in the stick-slip regime, tends to its maximum 
value F s tatic, and the relative variation in the stick-slip 
signal decreases. As a consequence, the role of thermally 
activated processes will be less strong for large Vq. More- 
over, the relative variation of Ff r i C with v s is much more 
pronounced for the lowest value of Vq, and the velocity 
dependence of friction becomes weaker for larger Vq. 




0.4 0.6 



1.8 2.0 



0.8 1.0 1.2 
v s [nm/s] 

FIG. 13. Friction force as function of scanning velocity for 
Vo = 2, Vo = 4 and Vo = 7. The filled circles connected by 
the dotted lines are the data for T = 0, while the open cir- 
cles connected by the dashed lines correspond to the data for 
T = 300 K. The solid lines are fits to the data at T = 300 K, 
according to Eq. (44) , with exponent a = 0.37 for Vo = 2, 
a = 0.56 for V = 4 and a = 0.67 ~ 2/3 for V = 7. The 
minimum value of the scanning velocity in the plot is v s = 1.5 
nm/s. 

The slope of the curves at T = 300 K slightly changes 
by increasing Vo and we find that the value 2/3 of the 
exponent of the logarithmic behavior (Eq. (42)) is recov- 
ered for the largest Vo we have used. This is in com- 
pliance with the approximation used to derive Eq. (42), 
namely Vo > 1 and Vo 3> k B T. More generally the data 
can be fitted by 



Ffric ~ F c oc | \n(v s 



(44) 



where the exponent a depends on Vo . In particular, from 
our data we obtain a(V = 2) = 0.37, a(V = 4) = 0.57 
and a(Vo = 7) = 0.67. A change of the slope of 
the velocity-friction curves can also be appreciated in 
Fig. 1(a) of Ref. [19], where data for different temper- 
atures are presented. This indicates that thermal effects 
critically depend on the surface corrugation and on tem- 
perature. 



V. DISCUSSION AND CONCLUSIONS 

In this paper we have investigated the velocity depen- 
dence of sliding friction at the atomic scale within the 
framework of the Tomlinson model. We have emphasized 
the role of the athermal processes characterizing the dy- 
namics, which are responsible for a power-law velocity 
dependence of the friction force at small scanning veloc- 
ities, while at finite temperature a creep regime takes 
place, giving rise to a logarithmic behavior of the friction 
force as a function of velocity. At variance with pre- 
vious claims in the literature, these dependences apply 
both in ID and 2D. We have also suggested in a semi- 
quantitative manner in which conditions thermal effects 
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are expected to be important for the frictional dynam- 
ics. Experimentally, the possibility to observe a veloc- 
ity dependent frictional force may crucially depend on 
the nature of the system, which determines the corruga- 
tion Vo, on the stiffness of the cantilever and on the ap- 
plied loading force, which in turns affects the value of V - 
Our model is simplified in the sense that the cantilever 
is treated as a point-like object and the form of energy 
dissipation, taken into account by introducing a damping 
term in the equations of motion, is purely phcnomcnolog- 
ical. Of course, in real situations finite contacts between 
the tip and the surface are involved and energy dissipa- 
tion comes into play through more complex mechanisms. 
However, a simple description based on the Tomlinson 
model contains the essential ingredients of the problem 
and can still capture the main dynamical features de- 
termining energy dissipation. We expect our study to 
stimulate further theoretical and experimental work on 
this issue. 
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